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Abstract 

The Karmen anomaly can be interpreted as being due to a heavy neu- 
trino of mass around 137 MeV and mean life 14 s produced in 7r e 2 decays. 
This interpretation is consistent with the present limits on the couplings 
of such an object. 
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1 What is the KARMEN anomaly ? 



The Karmen collaboration has reported a distortion in the time distribution 
of their events Q. The experiment studies low energy neutrinos produced by 
the dump of 800 MeV protons at the ISIS facility of the Rutherford Labora- 
tory. The accelarator has a unique time structure which allows to discriminate 
between neutrinos coming from n decaying at rest and from the subsequent \x 
decay at rest. It is in this last process that an anomaly was observed in the 
time distribution of interactions, which overall reproduces the lifetime of the 
muon but for an excess of events at a time delay centered on 3.6/xs w.r.t. the 
arrival of the protons. The signature is a deposition of energy in the range 11 
to 35 MeV. 

This is interpreted by the Karmen collaborators as possibly coming from the 
production of a heavy neutral object X of mass 33.9 MeV, travelling with a 
(3 = v/c of 1/60 and decaying in the detector, 17.5m away from the dump. The 
X particle is produced in tt decays : it — > [iX with a very small branching ratio 
and decays in turn into the a 'priori open channels: X — > or X — > ve + e~ or 
X — > z^77 

The Karmen experiment finds about 100 events compatible with X produc- 
tion and decay in a sample of 3000 neutrino interactions. This is a more than 
4 a effect. 

Experimentaly, there are two unknowns in the problem: the branching fraction 
Br of 7r's decaying into the X particle and the lifetime r of the latter. In order 
to reproduce the observed number of events, Br and r must be related by the 
curve shown in Fig.l as given in [|]|. For long lifetimes, the probability of decay 
in the detector is inversely proportionnal to the lifetime. For shorter lifetimes, a 
good fraction of X's has decayed before reaching the detector. Observe that, for 
very short r, the distribution of excess events in the detector should be uneven 
because of the stronger probability to decay close to the upper end. Taking this 
into account would change the shape of the curve in the low r region. However, 
Karmen sees no such effect, which can be taken as an indication that of the 
two possible r values corresponding to a given Br, only the higher one must be 
considered. 

A direct search for X production in ir^2 decay by looking at the muon spec- 
trum was performed at PSI [Q], ||. It excludes branching fractions larger than 
2.6 x 10~ 8 , as also indicated in Fig. 1. 
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2 Production and decay of massive neutrinos 



Massive neutrinos, arising, for example, from the introduction of right-handed 
chiral states, occur naturally in most extensions of the standard model. They 
have been searched through their decays in various experiments || 



2.1 Production 

If massive neutrinos exist, the so-called weak eigenstates (associated to charged 
leptons via weak currents) which we will call generically v\, are linear com- 
binations of mass eigenstates Vh ■ Let us call Uh\ the coefficient of Vh in the 
expansion of v\ (A being e or [i for the production processes envisionned here). 
In any V^' neutrino beam, there is a Vh component at the | ?1/X.A 1 2 level up to 
phase space factors and weak currents matrix elements. 

For example, 7r e 2 decays may yield Vh at the level \Uh e \ 2 up to the above 
mentionned factors, provided that has a mass below m(7r) — m(e) = 139 
MeV, while tt^ decays would produce Uh at the level \Uhfj,\ 2 with masses up to 
m(n) - m(fi) = 34 MeV. 

Note that the well known chirality suppression of 7r e 2 no longer occurs as soon 
as m v « a few MeV . 

For the X particle, A = fi and the mixing discussed in Fig 1 is . Since rax 
is supposed to be at the limit of the available phase-space in 71^2 (in order to 
ensure = 1/60), one would expect it to be produced more abundantly in K^. 
The ratio of the branching ratios ^~}^Z^^x) ~ ^ ins tead of .63 for massless 
v's. However, this is very sensitive to the precise value of rax- 



2.2 Decays 

Several decay modes can be investigated. Heavy neutrinos with masses below 
1 MeV can only decay radiatively by emitting one or two photons. As soon as 
the mass is above 1 MeV, the most favoured decay mode is e + e~v e . The matrix 
element is analogous to that of muon decay and the width is given by 

G 2 

T(v h -> e~e + v) = -^-^mlj (m Vh )\U he \ 2 

where / tends rapidly to 1 for masses above 10 MeV ([||]) and \Uh e \ is the same 
mixing which appears at production if is produced together with an electron. 
For higher masses, other decay channels open: ue^/j,^ 1 , e~ir + , h~tt + , ufi^fi + ... 
each one involving a mixing matrix element which depends on the final state 
leptons. 
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3 X or X* particle? 



The X particle could be such a massive neutrino undergoing weak decays. How- 
ever for a 33.9 MeV object e + e~v e is the only open non-radiative channel and 
its width, far above those of the radiative modes, dominates the lifetime. An 
upper limit of | t//,, e | 2 < 5 x 10~ 6 has been obtained Q for such a mass. This 
translates into a limit r > 120 s which is reported in Fig. 1 and excludes most 
of the allowed parameter region. For the highest mass eigenstate, published 
limits on r(vr -> ev)/r(ir -» \iv) (gj) allow to infer \U he \ 2 < 1.23 x 10" 6 and 
r > 500 s 

Moreover, taking the experimental Karmen curve (see Fig. 1) into account, this 
implies 

Br(-K + — > n + X) > 4.5 x 10 -9 which is barely compatible with the upper limits 
from the two PSI groups who tried a direct X search ||) [] 

An alternative interpretation has been proposed in terms of exotic muon 
decays. It involves unknown scalars. We study here yet another possibility, 
which implies an 'ordinary' massive neutrino, call it X*. 

As already remarked, many channels open as soon as phase-space permits. 
If we assume a 7r e 2 production mode, masses up to 139 MeV are possible, and 
ve^ur^ decays are allowed for sufficiently high mp, as schematically shown 
in Fig. 2 . For such modes, the time delay between beam on target and the 
events in excess in the Karmen detector would be the sum of the final state 
muon lifetime with the time necessary for X* to fly from the proton target to 
the detector. Indeed, the kinetic energy of /x^ is small and that of may fall 
short of the Karmen detection energy threshold, depending on nix* , rendering 
it necessary to wait for the muon to decay. 

This would yield an average (electron) energy around 36 MeV after an average 
delay of 2.2/xs following the undetected X* decay. 

The Karmen anomalous events appear to be rather tightly clustered in time if 
they are interpreted in the lines of Q and their apparent dispersion cannot be 
understood as due to a muon decay. However, the authors of || show that a 
rather narrow bump can be interpreted as the superposition of two exponen- 
tials, both with the same (muon) decay time constant, but with the second 
shifted by a more or less fixed amount representing a time of flight. 



x As quoted in the Review of Particle Physics Jj the PSI results would seem to contradict 
the lower bound we have just given and rule out altogether the Karmen solution. However 
the PDG are wrong in that they call limits on |C/x/j| 2 figures which are really limits on the 
branching ratio. Since, on the other hand, the way the PSI groups compute their upper 
bounds does not yield the tightest possible result, their is ample room for argument, and very 
little room, if any, for X 

2 for X* at rest, one finds E e < 14 MeV and T M < 2 MeV. 
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3.1 Model for the time distribution 



Therefore, we have tried to perform a two component fit to the Karmen time 
distribution. The first component is an exponential reflecting ordinary muon 
decay in the target followed by practically instantaneous interaction of the nor- 
mal light neutrinos thus produced in the detector (the ToF is about 60 ns) 
The second component is the convolution of a uniform and an exponential dis- 
tribution representing the sum of the transit time of a heavy neutrino from the 
target to the detector (which is also the heavy neutrino lifetime in the lab and 
is practically uniform in the window of interest) plus the lifetime of the muon 
produced in the heavy neutrino decay. 

However, this turned out to be inconsistent. The fitted ToF of the heavy 
neutrinos implies a mass close to 137 MeV which in turn means a positron 
energy distribution (in X* — > [i~e + v) rather evenly distributed around an ex- 
pectation value of 14 MeV, thereby invalidating the hypothesis that the delay 
would always comprise the \i lifetime. It was therefore necessary to add a third 
component representing those 'direct' positrons. Thus, there are altogether 5 
parameters representing the intensities of the three components and the lower 
and upper limits of the time window corresponding to the passage of the X* 7 s 
in the Karmen detector. 

In principle, this gives us a potential consistency check of the fact that the 
anomaly is really due to a particle since its speed can be estimated from the 
last two parameters using the average distance between the proton target and 
the detector on the one hand, and the length of the latter on the other hand. 
However, the data at our disposal are not accurate enough to decide between 
different fit results which do or do not conform with this expectation. Given 
the large errors on the fitted parameters, the results presented on Fig. 3 cannot 
be considered as inconsistent with the expectation. 

Here, Tl and T2 are the lower and upper edges of the time window, and the 
three following items are constants multiplying the three normalized distribu- 
tions superposed in the fit. Since the bin width is .5fis, the corresponding 
number of fj,~ — ► and X* — > e + are about 66 and 38 respectively. 
The fit yields: ToF = 3.66^s, hence (3 = .0159 and m x * ~ 137.3 MeV. 



3.2 Interpretation 

It must be remarked now that for an 'ordinary' heavy Dirac neutrino, there are 
really two (non mutually exclusive) possibilities: 

- either it decays to e" and a virtual W + which yields 

- or it decays to yT and the virtual boson goes into e + u e . 

The first possibility involves the same mixing matrix element, namely Uhe, 
which appears at production, whilst in the second case, the decay is due to 
Clearly enough, in any case Uhe 7^ and the e + e~v e mode is also widely open. 
If we assume the first channel, the unknown mixing cancels in the ratio of the 
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two modes which is given by the following formula: 



r = (1 _ 8r + r 2 )(1 _ r 2N _ i2rLog{r) 
T(e e+u e ) 

where r = (m^/mj.) 2 and m e = m Ve = m v = ||. The width ratio is about 
0.5% for m x * = 137 MeV. 

This shows two things: first, for the range of masses of interest, the lifetime 
is dominated by the e + e~v e mode, which gives us another relation between 
\Uhe\ and t and allows to draw further conclusions. Second, decays through 
this mode must have occured abundantly in the Karmen experiment, but the 
corresponding (large) deposition of energy is probably rejected by the Karmen 
selection. 

The results of this analysis are summarized on Fig. 4 It reproduces the Karmen 
'efficiency' curve corrected for the ir decay mode considered here (The unknown 
constant was taken from the asymptotic direction of the Karmen plot Q); also 
drawn is the branch of hyperbola corresponding to the Br x r value fixed by the 
X* mass and the hypothesis of an e~e + v dominated X* lifetime. The latter is 
found to be tx* ~ 14s and the ir + — > e + X* branching ratio is around 4 x 10 -10 
This entails a mixing matrix element |c7x*e| 2 ~ 4.2 x 10 -8 compatible with the 
previously published upper limits |4|, |5[ 

Here again, this object should be much more abundantly produced in kaon de- 
cays. The 'ratio of ratios' is around 100 in this case. 

Of course, we do not claim that this is "the" solution, but only wish to point out 
another possibility. It could probably be easily disproved through confrontation 
with more detailed Karmen data. 



4 Conclusion 



In short, the Karmen collaboration have chosen to explain their anomalous 
events by assuming they originate from a rare 71^2 decay, because a 'neutrino' 
with a mass at the edge of phase-space (around 34 MeV) has both a low enough 
velocity to explain the 3.6/is delay and a low enough mass to explain the 30 MeV 
energy deposition. We claim that both features can be explained by a 7r e 2 decay 
yielding a much heavier 'neutrino' which decays itself through a e^fj/^ui channel, 
the signature of which is a visible energy in the range of interest after the 
(apparently) observed delay. The more mundane e ± e^^ channel is supposedly 
hidden by the experimental energy cut off and the apparent clustering in time 
due to a fixed time of flight followed by a short detectability window is not 
spoiled upon taking the muon lifetime into account. 
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X* propagates from target to detector JL^ decays in detector 



Fig. 2 The process hypothetized in this paper 
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